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Abstract 

We observed an area of 10 deg 2 of the Large Magellanic Cloud using the Infrared Camera on board 
AKARI. The observations were carried out using five imaging filters (3, 7, 11, 15, and 24 micron) and 
a dispersion prism (2-5 micron, A/AA ~ 20) equipped in the IRC. This paper describes the outline of 
our survey project and presents some initial results using the imaging data that detected over 5.9xl0 5 
near-infrared and 6.4xl0 4 mid-infrared point sources. The 10 a detection limits of our survey are about 
16.5, 14.0, 12.3, 10.8, and 9.2 in Vega-magnitude at 3, 7, 11, 15, and 24 micron, respectively The 11 
and 15 micron data, which are unique to AKARI IRC, allow us to construct color-magnitude diagrams 
that are useful to identify stars with circumstellar dust. We found a new sequence in the color-magnitude 
diagram, which is attributed to red giants with luminosity fainter than that of the tip of the first red giant 
branch. We suggest that this sequence is likely to be related to the broad emission feature of aluminium 
oxide at 11.5 micron. The 11 and 15 micron data also indicate that the ([11] — [15]) micron color of both 
oxygen-rich and carbon-rich red giants once becomes blue and then turns red again in the course of their 
evolution, probably due to the change in the flux ratio of the silicate or silicon carbide emission feature at 
10 or 11.3 micron to the 15 micron flux. 

Key words: Galaxics:Magcllanic clouds infrared:stars stars:AGB and post-AGB, YSO ISM:supernova 
remnants 



1. Introduction 

Owing to its proximity (~ 50 kpc; e.g., Feast & Walker 
1987) and face-on geometry, the Large Magellanic Cloud 
(LMC) is an ideal natural laboratory for the study of var- 
ious astrophysical fields. It is located at a high Galactic 
latitude of ~ —36°, and we expect less contamination of 
foreground stars and less interstellar extinction. It is far 
enough to neglect its depth so that we can reasonably 
assume that objects in the LMC are all at the same dis- 
tance from us. The apparent size of the LMC is also in 
a size, for which the entire galaxy can be surveyed in a 
reasonable amount of time to study the material circu- 
lation processes and star-formation history in a galactic 
scale. In the meanwhile, it is close enough to resolve and 
study individual objects even with relatively small tele- 
scopes. Moreover, the mean metallicity of the LMC is 



known to be small (~ 1/4) compared to the solar, intrigu- 
ing us to study the influence of low metallicity on various 
astrophysical phenomena. 

Because of these characteristics there have been a num- 
ber of survey projects of the LMC in various wavelengths 
(Hard X-ray: Gotz et al. 2006; Soft X-ray: e.g., Long, 
Hclfand, & Grabelsky 1981 with Einstein, Sasaki et al. 
2000, Haberl & Petsch 1999 with ROSAT, and Haberl, 
Dennerl, & Pietsch 2003 with XMM; UV: Smith, Cornett, 
& Hill 1987; Ha: Kennicutt & Hodge 1986, Gaustad et 
al. 2001; Optical: Zaritsky et al. 2004; NIR: e.g., Cioni 
et al. 2000a with DENIS, Nikolaev & Weinberg 2000 
with 2MASS, and Kato et al. 2007 with IRSF/SIRIUS; 
MIR&FIR: e.g., Israel & Schwering 1986 with IRAS, Egan 
et al. 2003 with MSX, and Meixner et al. 2006 with 
Spitzer /SAGE; [CII]: Mochizuki et al. 1994; CO: e.g., 
Mizuno et al. 2001 with NANTEN; HI: Luks & Rohlfs 
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Fig. 1. Observed area of the AKARI IRC survey (thick solid 
outline), overlaid on the photographic image kindly provided 
by Mr. Motonori Kamiya. The thick dashed outline indicates 
the coverage of the Spitzer /SAGE survey (7° x7°, Meixner et 
al. 2006), the thick dash-doted outline shows the coverage of 
the IRSF/SIRIUS near-infrared survey (6.3°x6.3°, Kato et 
al. 2007), and the thin solid outline represents the coverage of 
the Magellanic clouds optical photometric survey (8.5° X 7.5°, 
Zaritsky et al. 2004). 

1992, Kim et al. 1998). 

A ground-based near-infrared (NIR) survey detected 
about fifteen million sources in the ~ 40 deg 2 area of the 
LMC (Kato et al. 2007). On the other hand, previous mid- 
infrared (MIR) surveys detected only a few thousands of 
sources in - 100 deg 2 area of the LMC (Egan ct al. 2003), 
which is too shallow to compare with the survey observa- 
tions in other wavelengths. Moreover, the angular reso- 
lutions of previous MIR surveys were not high enough to 
make secure cross-identifications with other survey cat- 
alogs. This situation is significantly improved with the 
advent of the Spitzer Space Telescope (SST; Werner et 
al. 2004) and AKARI, which have instruments capable of 
deep mid-infrared observations with the high spatial res- 
olution. The Spitzer SAGE project (Meixner et al. 2006) 
carried out an uniform and unbiased imaging survey of 
about 49 deg 2 area of the LMC, providing photometry for 
about 2xl0 5 sources with all IRAC (Fazio et al. 2004a) 
bands at two epochs separated by 3 months. 

We carried out near- to mid-infrared imaging and 
near-infrared spectroscopic survey toward the LMC with 
AKARI. The entire LMC was also observed as part of 
the All-Sky survey at 6 bands in the mid- to far-infrared 
with AKARI (Ishihara et al. 2006; Kawada et al. 2007). 
AKARI observations provide multi-band (11 bands) data 
of the LMC from the near- to far- infrared (FIR), which 
will give us a new insight into the various phenomena oc- 



curring in the LMC. In this paper, we provide an overview 
of our LMC survey project and present some initial results 
using the preliminary photometric catalog of bright point 
sources. 

Then we will make general analyses using the catalog 
with an emphasis on the unique characteristics of the IRC 
observations compared to Spitzer observations. The first 
AKARI/IRC LMC point source catalog is planned to be 
released to the public in 2009. 

2. The AKARI IRC LMC survey project 

The Japan Aerospace Exploration Agency launched an 
Infrared Satellite, ASTRO-F (Murakami et al. 2007) at 
21:28 UTC on February 21st, 2006 from the Uchinoura 
Space Center. Once in orbit ASTRO-F was renamed 
"AKARI". AKARI has a 68.5 cm telescope and two sci- 
entific instruments, namely the InfraRcd Camera (IRC; 
Onaka et al. 2007) and the Far-Infrared Surveyor (FIS; 
Kawada et al. 2007). Both instruments have low- to 
moderate-resolution spectroscopic capability. The IRC 
has nine imaging bands and six dispersion elements cov- 
ering from 2.5 to 26 /im wavelength range (Unfortunately, 
one of the dispersion elements that was expected to cover 
11 to 19 /im range became opaque during the ground test 
operation, and defunct in orbit). The FIS will observe 
in 4 far-infrared bands between 50 and 180 /im. One of 
the primary goals of the AKARI mission is to carry out 
an All-Sky survey with the FIS and the IRC at six bands 
from 9 to 180 /im (Ishihara et al. 2006; Kawada et al. 
2007). As a result, the entire LMC has been mapped in 
9, 18, 65, 90, 140, and 160 /im wavebands. The 5 o detec- 
tion limits for one scan observation of the All-Sky survey 
are shown in Figure 3. As described below, the LMC is 
located at a relatively high visibility region for AKARI's 
orbit. AKARI scanned LMC several times and thus bet- 
ter detection limits should be obtained by further data 
reduction of the All-Sky survey data. The All-Sky survey 
observations of the LMC will be described elsewhere and 
will not be discussed in this paper. 

In addition to the All-Sky survey in mid- and far- 
infrared, AKARI carried out two large-area legacy sur- 
veys (LS) in pointing mode. The LMC survey project 
(PI. T. Onaka) is one of the two LS programs. The other is 
the North Ecliptic Pole survey project (PI. H.Matsuhara; 
Matsuhara et al. 2007). These survey areas are located 
at high ecliptic latitudes, where the visibility is high for 
AKARI's sun-synchronous polar orbit. Therefore a batch 
of observing time can be allocated for pointing observa- 
tions in these areas after allocating scan paths for the 
All-Sky survey. We used the IRC to make imaging and 
spectroscopic mapping observations of the main part of 
the LMC. 

The IRC is comprised of three independent channels; 
NIR, MIRS, and MIRL, which cover the 1.8 - 5.5 /im, 
4.6 — 13.4 /xm, and 12.6 — 26.5 /im wavelength range, re- 
spectively. Each of the three channels has three imaging 
bands and two dispersion elements, which can be switched 
during a pointed observation opportunity. Each of the 
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Fig. 2. The normalized spectral response function of AKARI 
IRC bands and Spitzer IRAC and MIPS bands. As references, 
the ISO SWS spectra of two representative Galactic AGB 
stars (VX Sgr as O-rich AGB and IRC+50096 as C-rich AGB) 
with circumstellar dust features are shown. 

channels has a wide field of view (FOV) of about 10' x 
10', suitable to survey observations. Compared to the con- 
temporary SST the AKARI IRC has the following char- 
acteristics. 

• The IRC's imaging filters cover the wavelength 
range continuously from 2.5 to 26 /im. In particular, 
it has 11/xm (Sll) and 15/im (L15) imaging bands, 
which fill the gap between the IRAC (Fazio et al. 
2004a) and MIPS (Rieke et al. 2004) on board the 
SST. The wavelength range that Sll and L15 bands 
cover contains interesting spectral features such as 
the silicate 10 and 18 ^m bands. Figure 2 shows the 
spectral response functions of the IRC, IRAC, and 
MIPS bands, together with ISO-SWS spectra (Sloan 
et al. 2003) of an oxygen-rich AGB star (VX Sgr) 
and a carbon-rich AGB star (IRC +50096). The 
IRC Sll and L15 probe the silicate bands efficiently. 

• One of the most significant features of the IRC is its 
ability to perform slit-less spectroscopy in addition 
to imaging. It can simultaneously obtain near- to 
mid-infrared (about 2 — 13 /im) continuous spectra 
of all sources present in its large FOV of about 10' 
x 10' within a pointing opportunity. 

These features make AKARI IRC unique and complemen- 
tary to IRAC and MIPS on the SST. Refer to Onaka et al. 
(2007) for instrumental details and imaging performance 
of the IRC. General information on the IRC spectroscopic 
mode, particularly on the slit-less spectroscopy, is given 
in Ohyama et al. (2007). 



3. Scientific objectives 

Interstellar matter in a galaxy is thought to evolve 
through various processes. Stars are born in interstellar 
clouds by accreting the interstellar matter and then ap- 
pear as young stellar objects (YSOs), such as T Tauri stars 
or Herbig Ae/Be stars. Low- and intermediate-mass stars 
evolve into white dwarfs after spending most of their lives 
in the main sequence. During the course of their evolution 
to white dwarfs, they lose mass through so called mass-loss 
process and supply dust grains and gas to the interstellar 
space. High-mass stars end their lives as supernovae and 
enrich the interstellar medium with a large amount of nu- 
clear synthesized elements. They can also trigger the for- 
mation of next generation stars through their strong UV 
photons and stellar winds (Hosokawa & Inutsuka 2006; 
Koo et al. 2008). Interstellar grains play an important role 
in every aspect of these interstellar processes, including 
the star-formation. All of these phenomena can be stud- 
ied most effectively in the infrared. The Spitzer SAGE 
project (Meixner et al. 2006) have provided IRAC (3.5, 
4.5, 5.8, and 8.0 /Ltm) and MIPS (24, 70, and 160 /imj 
images for a wide area (7° x 7°) of the LMC. Adding to 
the SAGE observation, the complementary AKARI IRC 
(3.2, 7.0, 9.0, 11.0, 15.0, and 18.0 /im) and FIS (65, 90, 
140, and 160 /im) observations of the LMC will provide 
an opportunity to make a thorough study of these mate- 
rial circulation processes and local star-formation history 
of the LMC in a galactic scale owing to its wide spectral 
coverage and high sensitivities together with the sufficient 
spatial resolution. 

Figure 3 is a diagrammatic representation of 5 a detec- 
tion limits of our survey, Spitzer SAGE survey, and the 
IRSF/SIRIUS near- infrared survey. The 5 a detection 
limits of our survey are calculated by scaling the 10 a de- 
tection limits in Table 1 (see section 5.4.1). We take the 5 
a detection limits for the SAGE survey from the " SAGE 
Data Description: Delivery l" 1 . The detection limits are 
expected to be further improved in the final catalog deliv- 
eries (Meixner et al. 2006). It is clear that the detection 
limits of our survey are comparable to that of the Spitzer 
SAGE survey. The ISO SWS spectra of galactic famous 
stars of known distances (the distances were taken from 
Crosas & Menten 1997 for IRC +10216, and Perryman et 
al. 1997 for others) are included for illustrative purposes 
after scaling their fluxes at the distance of the LMC. We 
assumed a distance modulus of 18.5 mag to the LMC. 
The figure indicates that all of red giants above the tip of 
the first red giant branch, and a large fraction of Herbig 
Ae/Be stars in the LMC are detected in our survey. 

3.1. Star-formation and young stars 

The LMC contains a class of so called populous clus- 
ters (e.g., Kumai, Basu, & Fujimoto 1993) that have 
no Galactic counterpart. They have masses intermedi- 
ate between Galactic globular clusters and open clusters. 
Characteristic environments in the LMC, such as low 

1 http:/ /sage. stsci.edu/index.php 
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Fig. 3. A graphic representation of the 5 a detection limits 
of the AKARI LMC survey and the All-Sky survey (magenta). 
For comparison, the 5 a detection limits of the Spitzer SAGE 
survey and the IRSF/SIRIUS survey are also shown (cyan 
and gray, respectively). Some ISO SWS spectra (Sloan et al. 
2003) of famous galactic stars are scaled to fit the distance to 
the LMC and included in the figure, as examples of Herbig 
Ae/Bc stars (AB Aur), red giants with luminosities below the 
tip of the first red giant branch (/? And), and dusty red giants 
(IRC +10216) in the LMC. 

mean metallicity (about 1/4 of solar vicinity), low dust 
to gas ratio (about 1/4 of Milky Way; Koornneef 1982) 
and strong ultraviolet interstellar radiation (e.g., Cox & 
Spaans 2006), may play a key role in the star-formation 
process of populous clusters (Fukui 2005). Detailed stud- 
ies of the star-formation in the LMC will provide insights 
into these processes. The IRAS survey has detected pro- 
tostars in the LMC only brighter than 10 4 ~ 10 5 L Q (e.g., 
van Loon et al. 2005). Although only the brightest end 
of the YSOs in the LMC can be detected by the IRC (see 
the detection limits provided in Table 1), the IRC imaging 
data provide us with a valuable dataset to investigate the 
nature of such objects from the spectral energy distribu- 
tion from near- to mid-infrared as well as the water ice 
feature at 3 /an and/or silicate absorption at 10 and 20 
fxm. AKARI observations of the LMC enable us to probe 
the luminosity function of protostars down to about 10 
L0. 

The detection limits of IRC bands are not deep enough 
to detect T Tauri stars in the LMC. However the IRC 
observations can nearly completely detect young stars of 
intermediate mass (Herbig Ac/Be) as well as classical Be 
stars. It is difficult to separate Herbig Ae/Be stars from 
classical Be stars and another kind of young stars with dif- 
ferent masses and ages, based on near-infrared data alone. 
Mid-infrared data provided by the IRC make clear iden- 
tification of these objects. AKARI IRC observations are 
able to provide the initial mass function of objects later 
than the Class 1 phase for M > 2 - 3 M except for 
dense clusters, where the source confusion limits the de- 
tection. They allow us to compare the spatial distribu- 



tion of YSOs with that of CO clouds in a statistical way 
and investigate the star- formation history in the LMC. 
Particularly interesting is the search for embedded star 
clusters in " cluster less" giant molecular clouds (GMCs), 
such as N159S (Bolatto et al. 2000) and investigation of 
whether or not these clusterless GMCs are really not as- 
sociated with star clusters. The 5 band IRC observations 
together with 4 band FIS All-Sky survey data allow us 
to distinguish these YSOs unambiguously and study the 
nature of detected objects in detail. 

3.2. Evolved stars 

Stars of low- to intermediate-mass (stars with main 
sequence masses M ms of about 0.8-8 M Q ) eventually 
evolve into white dwarfs. The mass M W( j distribution of 
white dwarfs peaks around 0.6 M Q (e.g., Kepler et al. 
2007). Therefore, the difference, AM = (M ms - M wd ) 
should be lost during their life in the mass-loss process, 
which enriches gas and dust in the interstellar medium. It 
is observationally known that the maximum mass-loss rate 
is attained at the latest stage of their Asymptotic Giant 
Branch (AGB) evolutionary phase, where stellar pulsation 
is usually associated. However, the onset and evolution of 
the mass-loss process and its relation to pulsation are not 
yet completely understood. 

As described above, the LMC has a number of young 
clusters, which also provide a unique place to study the 
late evolutional stage of intermediate-mass stars. Star 
clusters of intermediate ages still contain a number of 
AGB stars that are losing mass. Since their ages and 
masses are well determined, the study of mass-losing AGB 
stars in the clusters can elucidate the mass-loss process for 
stars with given ages and masses quantitatively for the 
first time. 

Blum et al. (2006) used Spitzer SAGE survey epoch 
1 catalog to discuss and identify the evolved stars 
in the infrared color-magnitude diagrams. Ita et al. 
(2004a, b) studied OGLE variables in the Magellanic 
Clouds (Udalski, Kubiak & Szymanski 1997; Zebruh et al. 
2001). They detected most of optical variable stars in the 
central parts of the Large and Small Magellanic Clouds. 
The association of the IRC data with these variability cat- 
alogs enables us to investigate the mass-loss phenomena 
in terms of the pulsation activity directly. The IRC obser- 
vations can detect all mass-losing AGB stars in the LMC. 
It should be emphasized that the IRC band Sll, L15, and 
L24 are particularly important to determine the chemistry 
of the envelope (O-rich or C-rich) and the mass-loss rate 
accurately because of the presence of the silicate bands at 
10 and 18 /im and the silicon carbide band at 11.3 /im. 

3.3. Supernova remnants 

Supcrnovac (SNe) are believed to be a major source of 
interstellar dust since the expected supply rate from late- 
type stars is too short to balance with the dust destruction 
rate in the interstellar space (Jones et al. 1994; 1996). The 
role of the dust production in SNe is crucial to understand 
the circulation and evolution of materials in the inter- 
stellar medium (Dwek 1998). However, to date, no clear 
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observational evidence has been given for the large pro- 
duction of dust grains from SNe (Meikle et al. 2007). The 
condensation efficiency could be cither very low (Clayton 
et al. 2001) or moderate (Todini & Ferrara 2001). A study 
of Cas A SNR by Spitzer clearly demonstrates the obser- 
vational difficulty of dust grains associated with SNRs in 
our Galaxy (Krause et al. 2004). SNRs in the Galaxy are 
easily confused by the material on the line of sight because 
most of them are located on the Galactic plane. In this 
regard, the LMC provides an ideal opportunity to study 
material ejected in SNRs without suffering from serious 
confusion problems. Recent studies of SNRs in the LMC 
with Spitzer and AKARI provide significant information 
on the origin of the IR emission in SNRs and the inter- 
stellar dust lifecycle in terms of dust formation by SNe 
and dust distruction by SNR blast waves (e.g., Borkowski 
et al. 2006; B. Williams et al. 2006; R. Williams, Chu, & 
Gruendl 2006; Rho et al. 2008; Seok et al. 2008). 

The interaction of SNRs with their surroundings is 
also an important phenomenon, providing opportunities 
for studying strong shocks and sometimes triggered star- 
formation (Koo et al. 2007; 2008). However, the ob- 
servational study of the interaction is also hampered by 
confusion for Galactic SNRs. There are more than 40 
SNRs known in the LMC (Filipovic et al. 1998; Haberl & 
Petsch 1999), and only 10 of them have been detected by 
IRAS. SNRs appear very patchy in the infrared (Douvion 
et al. 2001) and their spectral energy distribution indi- 
cates the multi-temperature nature. They also have many 
emission lines, which could affect the broad-band photom- 
etry. Strong emission lines expected in the MIR are [Aril] 
7.0 /mr, [Nell] 12.8 /mi, [Nelll] 15.6 /im, [Fell] 17.9 fj,m, 
[Fell] 24.5, 26.0 fim, and [Sill] 34.8 /jm (e.g., Arendt et al. 
1999; Tcmim et al. 2006; Rho et al. 2008). Dust grains in 
SNRs may be transiently heated by high-energy electron 
hits, which results in a wide dust temperature distribu- 
tion (Dwek 1986). A peculiar band may also be present 
in the MIR (Arendt et al. 1999), which might be associ- 
ated with dust grains formed in the SNR (Chan & Onaka. 
2000; Onaka et al. 2008). Therefore, multi-band photom- 
etry spanning over a wide spectral range with sufficient 
spatial resolution is crucial to study the nature of the in- 
frared emission from SNRs and could confirm the physical 
association of the emission with SNRs. 

3.4- Interstellar Matter 

Diffuse infrared emission from interstellar dust consists 
of several components: thermal emission in the FIR from 
submicron size grains, excess continuum emission in 15.80 
/im, and the unidentified infrared (UIR) emission bands 
arising from aromatic materials (Onaka 2000). Each com- 
ponent should play a significant role in the interstellar pro- 
cess and their relative variations in abundance should pro- 
vide valuable information on the physical conditions of the 
interstellar medium. However, little has been known for 
formation, evolution (processing), and destruction of each 
dust component. A study of the UIR bands in our Galaxy 
has suggested possible variations in the relative abundance 
to the submicron grains as well as in the band appearance 



with the Galacto-centric distance for the first time (Sakon 
et al. 2004). Observations of the Galactic plane pick up 
diffuse emissions from every component on the line of sight 
and it is impossible to separate out the emission from a 
particular region. The LMC again provides an ideal place 
to study the diffuse infrared emission from dust grains 
in several different environments compared to our Galaxy 
because of its nearly face-on geometry. It has been sug- 
gested that the nature of the infrared emission from the 
LMC is different from our Galaxy (Aguirre et al. 2003), 
which may be related to the low metallicity of the LMC. 
A recent study of IRAS data has indicated that the color 
of the IRAS 12 /im to 25 /zm intensities seems to be high 
in CO clouds without young star clusters, whereas it is 
low in regions with YSOs (Sakon et al. 2005). There is 
also an indication that the 12 /im to 100 /jm ratio is differ- 
ent inside the supcrgiant shells (Kim ct al. 1998), which 
is also related to the star-formation activity (Yamaguchi 
et al. 2001). Since the IRAS 12 /im band intensity in- 
cludes both the UIR and the MIR excess emissions, it is 
important to separate the contribution from each compo- 
nent and compare it directly to the FIR component. The 
IRC S7 filter is most sensitive to the 6.2 and 7.7 /jm UIR 
bands, whereas the Sll band will probe the UIR 11.2 /im 
band. The S7 to Sll band ratio together with the MIR 
and FIR All-Sky survey data can investigate variations in 
the relative band strength efficiently The L15 and L24 
band data determine the contribution and spectral pro- 
file of the MIR excess component unaffected by the UIR 
bands. 

4. Observations 

4--1- IRC astronomical observing template and pointing 
observation seguence 

To cover a wide range of the spectral energy distribu- 
tion of a celestial body, we use the AKARI IRC AOT02 
observing template with a special option prepared only 
for the LMC survey. The AOT02 with the special option 
yields not only imaging data at 3 (N3), 7 (S7), 11 (Sll), 
15 (L15), and 24 /im (L24), but also low resolution (A/AA 
~ 20) 2.5 — 5 /im (NP) spectral data at three dithered 
sky positions in a pointing opportunity. Figure 4 shows an 
example of the NP spectroscopic data taken in our survey. 
The reduced spectrum of a red source indicated in the box 
is also shown in the figure. Discussions on spectroscopic 
data will be given in a future paper. 

As described in the IRC data user's manual (Lorcnte 
et al. 2007) the IRC observation produces short- and 
long-exposure data together. The long to short exposure 
time ratios in each channel are 9.5, 28, and 28 for NIR, 
MIRS, and MIRL, respectively. It enhances the dynamic 
range by combining photometric results of the two differ- 
ent exposure-time data. In this paper, we use the com- 
bined data of long- and short-exposure for N3, but only 
use long-exposure data for MIRS and MIRL. 
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Fig. 4. A sample image of NP slit-less spectroscopy data 
taken in our survey and a reduced NP spectrum of a source 
indicated in the box. 

4-2. Mapping strategy 

Due to the focal-plane layout of AKARI, the NIR and 
MIRS channels share the same field of view (i.e., they ob- 
serve the same sky position), but the MIRL channel ob- 
serves a different sky position that is offset by 25 arcmin 
from the N1R/MIRS center in the direction perpendicu- 
lar to the AKARl's orbit. AKARI is in a sun-synchronous 
polar orbit at 700 km attitude along the twilight zone, and 
its orbit is approximately parallel to the ecliptic meridian 
lines with approximately 4.f arcmin spacing at the eclip- 
tic plane between successive orbits. We used this feature 
to map the LMC effectively. Each of the IRC channels 
has a field of view of f0' x 10' and adjacent fields are 
designed to overlap by ~ 1 f 5 . Since AKARI follows the 
Earth's yearly round, the telescope points at the same ce- 
lestial position every six months, but the direction of the 
satellite movement relative to the position is rotated by 
180 degrees. Thus the position of the MIRL channel rel- 



ative to the NIR and MIRS channels is also rotated by 
180 degrees in projection on the sky. Therefore, we divide 
the observing area into several parts to maximize the ob- 
serving efficiency to map a large part of the LMC with 
all three channels of the IRC. Observations were carried 
out in three separate seasons, from 6th May 2006 to 8th 
June 2006, from 2nd October 2006 to 31st December 2006, 
and from 24th March 2007 to 2nd July 2007. Over 600 
pointing observations were devoted for this project, yield- 
ing about a 10 deg 2 imaging and spectroscopic map of the 
main part of the LMC. 

5. Data reduction 

5.1. Pipeline processing 

Raw imaging data were processed with the IRC imag- 
ing toolkit, version 20071017 (see IRC Data User's Manual 
Lorente et al. 2007 for details). In a pointing observation 
using IRC, one short-exposure and three long-exposure 
dark data are taken at the beginning (prc-dark) and the 
end of the operation for all of the three IRC channels. 
We made long-exposure dark images for MIRS and MIRL 
channels by averaging over three long exposure dark data. 
We refer to these averaged pre-dark images as the self- 
dark. We subtracted selfdarks from MIRS and MIRL 
long-exposure images and the pipeline default super-dark 
from the others. 

As described above, we use the AKARI IRC AOT02 
observing template that produces data at three different 
sky positions in a pointed observation. The data in each 
band are spatially aligned and then coadded by taking 
their median to eliminate array anomalies such as bad, 
dead, or hot pixels, cosmic rays, and so on. The resultant 
coadded images have pixel sizes of I. "446 pixel -1 , 2 "340 
pixel" 1 and 2 "384 pixel" 1 for NIR, MIRS, and MIRL im- 
ages, respectively. 

5.2. Astrometry 

After the coaddition process, we calculate the coordi- 
nate transform matrix that relates the image pixel coordi- 
nates to the sky coordinates by matching detected point 
sources with the Two Micron All Sky Survey (2MASS; 
Skrutskic et al. 2006) catalog. We use at least five 
matched point sources for the calculation. If matching 
with the 2MASS catalog is unsuccessful (such cases usu- 
ally occur in LI5 and L24 images), then we use the SAGE 
point source catalog (Meixner et al. 2006) as the positional 
reference. The root-mean-squares of the residuals between 
the input 2MASS/SAGE catalog coordinates and the fit- 
ted coordinates are smaller than 1"2, 2."6, and 2"9, for 
NIR, MIRS, and MIRL images, respectively. The coordi- 
nates of AKARI sources should be accurate to that extent 
relative to the 2MASS and SAGE catalog coordinates. 

5. 3. Photometry 

Photometry is carried out for each coadded image in- 
dependently. To derive calibrated fluxes for each detected 
source, point spread function (PSF) fitting photometry 
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was performed on the coadded images with the IRAF 2 
package DAOPHOT. We develop PSF fitting photome- 
try software, which is similar to the one that the Spitzer 
GLIMPSE 3 (Benjamin et al. 2003) team uses. Our pho- 
tometric process involves the following steps: 

1. Source extractor (Bertin & Arnouts 1996) is used for 
each coadded image to extract sources whose fluxes 
are at least 3 a above the background. The sat- 
urated sources should not be extracted, since the 
central pixels of them are already masked by the 
IRC imaging toolkit. Even if they are detected by 
the source extractor, the PSF fitting does not work 
well for them, and they are not included in the fol- 
lowing discussions. 

2. Aperture photometry is performed on all of the 
sources found in step 1, using the task PHOT with 
aperture radii of 10.0 and 7.5 pixels for NIR and 
MIRS/MIRL images, respectively. We use the same 
aperture radii as have been used in the standard star 
flux calibration (Tanabe et al. in preparation), and 
the aperture corrections are not applied. The resul- 
tant arbitrary fluxes are converted into the physical 
units by using the IRC flux calibration constants 
version 20071112. Then the calibrated fluxes are 
converted into the magnitudes of the IRC- Vega sys- 
tem by using the zero magnitude fluxes listed in 
Table 1. The offsets between the arbitrary and the 
calibrated magnitudes are constants. 

3. Several point sources with moderate flux (i.e., with 
a good signal-to-noise ratio and unsaturated) and 
without other sources within 7 pixels are selected 
from the results of step 2. At least 8 such "good" 
stars are selected in the NIR/MIRS (i.e., N3, N4, 
S7, and Sll) images, and more than 5 good stars 
in the MIRL (i.e., L15 and L24) images. They are 
used for the aperture correction to the PSF fitting 
photometry in the following process. 

4. Since the shapes of the PSFs measurably vary from 
pointing to pointing in NIR images possibly due to 
jitters in the satellite pointing, the selected good 
stars in step 3 are used to construct a model PSF 
for each image. We use DAOPHOT to choose the 
best fitting function by trying several different types 
of fitting functions. In MIRS/MIRL images, we 
use a "grand-PSF" for each band, which was built 
in advance. These grand-PSFs were made from 
"good" stars free from diffuse emission, carefully 
chosen by eyes from hundreds of images taken in 
the AKARI/IRC LMC survey data and other re- 
lated ones. 

5. The PSF fitting photometry is performed on all 
of the sources found in the coadded images using 

2 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 

3 Visit http://www.astro.wisc.edu/glimpsc/, and see a docu- 
ment "Description of Point Source Photometry Steps Used by 
GLIMPSE", written by Dr. Brian L. Babler. 



ALLSTAR. The PSF does not vary appreciably over 
the field of view. Hence a constant PSF is assumed 
over an image. This PSF fitting operation is iter- 
ated, in just the same way as the GLIMPSE team 
does. 

6. Aperture correction is applied to the fit magnitudes 
by comparing them with the aperture magnitudes of 
good stars, which are selected in step 3. Then, the 
offset value obtained in step 2 is added to the aper- 
ture corrected fit magnitudes to derive calibrated fit 
magnitudes. 

As indicated above, we use the zero magnitude fluxes 
(Tanabe et al. in preparation) to convert the photometric 
fluxes into the Vega magnitudes. These zero magnitude 
fluxes are tabulated in Table 1. It should be noted that the 
IRC absolute flux calibration assumes a spectral energy 
distribution (SED) of fx oc A -1 or f v oc v~ x . Note that 
throughout this paper, we do not apply color corrections 
to the calibrated fluxes and magnitudes. No dereddening 
are applied either. 

5-4- Compilation of a preliminary catalog 

We carried out photometry on each coadded image as 
outlined in the previous section. It produced photometric 
catalogs for each band for individual images. At first, the 
catalogs for individual IRC imaging bands (i.e., N3, S7, 
Sll, L15, and L24) are made. Because of the field over- 
lap, there should be substantial multiple entries of a single 
source in the catalogs. We eliminate such multiple entries 
based on the spatial proximity (| Ar| < 1"5). We adopt the 
one with a better S/N and discarded the other(s). We use 
both long- and short-exposure data of N3, and we merged 
them using a positional tolerance of l."5. If a source is 
detected in both short and long exposure data, we adopt 
the one with a better S/N and discard the other. Then, 
two corresponding band catalogs are merged to make in- 
dividual IRC channel catalogs. The, S7 and Sll (L15 and 
L24) catalogs are merged using a positional tolerance of 
l."5 (l."5) to make a MIRS (MIRL) channel catalog. For 
the matched sources, their coordinates are recalculated by 
taking an average of the coordinates of each band. In rare 
cases that more than two sources are present within the 
tolerance radius, the closest one is always adopted and 
the others were listed as solitary sources. The resultant 
individual channel catalogs are further merged to make a 
grand catalog using a positional tolerance of 1"5. 
5.4-1- Detection limits and completeness 

The distributions of the photometric uncertainty versus 
the magnitude in the IRC bands are shown in Figure 5. 
The photometric uncertainty includes errors in the fitted 
magnitude calculated by IRAF/ ALLSTAR, in the aper- 
ture correction factor, and also in the ADU-to-Jy conver- 
sion factor. The horizontal dashed lines show the signal 
to noise ratio (S/N) of 10, and the vertical solid lines show 
the 10 a detection limits, which are defined as the faintest 
magnitudes at which the mode of photometric uncertain- 
ties of the sources in 0.02 mag bins exceed 10 a, except for 
the N3 data. In the calculation of the 10 a detection limit 
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1 These are the most common values. Total exposure time per pixel can be different from coadded images to images, 
because some data badly damaged by cosmic-ray and/or hitting of charged particles are discarded. 

2 For point sources. 

3 Numbers arc taken from ASTRO-F Observer's Manual version 3.2 (ASTRO-F Observer's Manual. 2005). 

4 Tanabe et al. in preparation. 



for the N3 data, we use only data detected in long expo- 
sure images, and take a mean of photometric uncertainties 
in 0.02 mag bins instead of taking a mode of them. Visual 
inspections showed that sources can be detected only in 
short exposure images, and this is why there is a sequence 
around 15 mag in the N3 error plot. There seems to be 
several causes to explain this. For example, not a small 
number of long exposure images are damaged by column 
pulldowns caused by saturated stars, and sources in such 
polluted columns are not detected in the long exposure 
images. The estimated 10 a detection limits are 16.50, 
14.00, 12.54, 10.74, and 9.16 in N3, S7, Sll, L15, and 
L24, respectively. These detection limits are comparable 
to that of the Spitzer SAGE survey (see Figure 3). 

The magnitude distributions of the sources in the pre- 
liminary catalog are shown in Figure 6. The sizes of the 
bins are 0.1 mag for N3 and 0.2 mag for S7, Sll, L15, 
and L24. The vertical solid lines show peak magnitudes, 
below which photometry is incomplete. The peak magni- 
tudes are 15.9, 12.8, 12.2, 10.4, and 8.8 mag in N3, S7, 
Sll, L15, and L24, respectively. 

The 10 cr detection limits for point sources and source 
counts in each IRC band are summarized in Table 1 to- 
gether with other properties of our survey. 
5.4-2. Cross-identifications with the existing point source 
catalogs 

We cross-identified our catalog with the following exist- 
ing catalogs using a positional tolerance of 1"5, and use 
the result for discussion in the rest of the paper. 

• The Magellanic Clouds Photometric Survey catalog 
(Zaritsky et al. 2004): It lists U, B, V, and I stellar 
photometry of the central 64 deg 2 area of the LMC. 

• The IRSF Magellanic Clouds Point Source Catalog 



(Kato et al. 2007): The IRSF catalog lists JHK S 
photometry of over 1.4xl0 7 sources in the central 
40 deg 2 area of the LMC. Compared to the con- 
temporary DENIS (Cioni et al. 2000a) and 2MASS 
(Skrutskie et al. 2006) catalogs, the IRSF catalog is 
more than two mag deeper at K s -band and about 
four times finer in spatial resolution. 
• The Spitzer SAGE catalog (Meixner et al. 2006): It 
lists near- to far-infrared photometry of sources in 
the central 49 deg 2 area of the LMC. 

The observed regions of the above existing catalogs are 
indicated in Figure 1. The figure shows that a part of the 
AKARI IRC survey region is outside the bounds of the 
IRSF/SIRIUS survey area. We do not use the AKARI 
sources located outside the IRSF/SIRIUS area in the fol- 
lowing discussions. Throughout this paper, the numbers 
bracketed by [ ] designates the data of the SAGE catalog, 
for example, [3.6] indicates the photometry in the IRAC 
3.6 /zm band. 

5-4-3. Comparison to the Spitzer SAGE survey catalog 

Although the bandpasses of the AKARI IRC and the 
Spitzer IRAC/MIPS bands are different, comparison of 
closely matched bands is useful to test the calibration of 
the AKARI IRC data. We compared the IRC N3, S7, 
L24 photometries in our catalog with the corresponding 
3.6, 8.0, and 24 //m fluxes of the sources in the Spitzer 
SAGE catalog. Note that both of the IRC and IRAC 
absolute flux calibration assume a SED of /„ oc v" 1 , but 
MIPS flux scale assumes a source spectrum of a 10,000 
K blackbody (MIPS Data Handbook 2007). Therefore 
we converted MIPS scale into IRC/IRAC one by adding 
-0.043 mag (Bolatto et al. 2006) to the MIPS 24 cata- 
log magnitudes to make a direct comparison between L24 
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Fig. 5. Photometric uncertainties as a function of the mag- 
nitude at each IRC band. The sizes of the bins are 0.1 X 
0.01 mag for N3, and 0.2 X 0.01 mag for S7, Sll, L15, and 
L24. The dashed lines show 10 a errors, and the solid lines 
indicates the 10 a detection limits: 16.50, 14.00, 12.54, 10.74, 
and 9.16 mag in N3, S7, Sll, L15, and L24, respectively. 

and MIPS 24 photometries. The distributions of the 
magnitude differences between the IRC and IRAC/MIPS 
fluxes as a function of the corresponding IRC magnitudes 
are shown in Figure 7. The triangles show the mean of 
the residual magnitudes in the corresponding 1 mag bin, 
and the error bars show their 1 a standard deviations. 
Taking account of the differences in the bandpass, the 
AKARI IRC and IRAC/MIPS photometries appear con- 
sistent with each other. Quantitatively, they are in agree- 
ment within 9, 10, and 11 % in N3 vs [3.6], S7 vs [8.0], and 
L24 vs [24] comparisons, respectively, using the sources 
with S/N > 10 a and after subtracting systematic offsets. 

In the comparison between N3 and [3.6] photometries, 
we find an interesting feature. It is seen that a group 
of bright sources (sources with 9 < N3 [mag] < 11) devi- 
ate from the overall trend. To explain the deviation, we 
cross-identified our sources with spectroscopically identi- 
fied carbon stars (Kontizas et al. 2001), which are shown 
in blue dots. We photometrically select dusty carbon-rich 
AGB candidates that cannot be identified in the optical 
survey based on the following conditions. 




6 8 10 12 4 6 8 

Apparent Magnitude [mag] 

Fig. 6. Magnitude distribution of the sources in the AKARI 
LMC survey catalog. The sizes of the bins are 0.1 mag for 
N3, and 0.2 mag for S7, Sll, L15, and L24. The solid lines 
show the peak of the source count histograms, below which 
the photometry is incomplete: 15.9, 12.8, 12.2, 10.4, and 8.8 
mag in N3, S7, Sll, L15, and L24, respectively. 



• (J - K s ) > 2.0 

• K s < (4/3) x ((J 



K s ) - 2.0) + 12.0 



These sources are indicated by red dots. The J and K s 
fluxes are obtained by cross-identifying the IRC catalog 
with the IRSF/SIRIUS LMC near- infrared catalog (Kato 
et al. 2007, see next section). It is true that the above 
conditions may include some dusty O-rich AGB stars, but 
such stars are small in number compared to the carbon- 
rich AGB stars. It is clearly seen that red dots are the 
sources that show the deviation in the diagram. The de- 
viation should be attributed to the fact that the passband 
of the IRC N3 band is slightly bluer than that of the IRAC 
3.6 ^m, and it is more sensitive to the 3.1 /mi HCN-I-C2H2 
absorption feature, which is usually seen in carbon-rich 
AGB stars (e.g., Yamamura & de Jong 2000). The N3 
data, thus, enable us to select sources with carbon-rich 
circumstcllar chemistry on the basis of their colors. 
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Table 2. Source Classification based on Kraemer et al. 2002. 
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Fig. 7. Comparison between AKARI LMC survey and 
Spitzer SAGE catalog using closely matched passbands (top 
panel: N3 vs [3.6], middle panel: S7 vs [8.0], bottom panel: 
L24 vs [24]). In the N3 vs [3.6] panel, the spectroscopically 
identified carbon stars (optical prism survey, Kontizas et al. 
2001) are indicated by blue dots, and photometrically selected 
(see text for selection conditions) carbon-rich AGB candidates 
arc shown in red dots. 



6. General analysis using the AKARI LMC sur- 
vey photometric catalog 

6.1. Color-color diagram 

6.1.1. Convolution of ISO SWS spectra 

With an aim to help an interpretation of what types 
of sources are detected by the AKARI survey, we take 
fully calibrated ISO SWS spectra data from the on-line 
database 4 available (Sloan et al. 2003). Then we estimate 
their fluxes in the AKARI IRC system by convolving the 
ISO SWS spectra with the response curves of the IRC 
bands using the following equation, 



/•reference 



X > hv 



(1) 



where the suffix i denotes IRC bands, A.; represents the 
reference wavelengths of each IRC band tabulated in ta- 
ble 1, and Ri(X) indicates the spectral responses of each 
band. The calculated fluxes are further converted into the 
Vega-magnitudes. Note that the calculated magnitudes 
are in the IRC system (substituting A/a (A) = constant to 

equation (1) gives Ai/^ e ^ erence ), which can be directly 
compared with the catalog values. 

We divide ISO SWS sources into 10 main groups, based 
on the Table 6 in Kraemer et al. (2002). Sources that are 
classified as "uncertain classification" and "very uncertain 
classification" in Kraemer et al. (2002) are excluded in the 
following discussions. Table 2 summarizes how we group 
the ISO SWS sources. 

6.1.2. (Sll - L15) vs (N3 - Sll) diagram 

Figure 8 shows the color-color diagrams of AKARI 
sources (upper panel) and ISO SWS sources (lower panel) 
using colors of (Sll - L15) and (N3 - Sll). In the 
AKARI diagram, only stars with S/N > 5 in all employed 
three bands are included. The numbers of sources in each 
0.05 x 0.05 mag 2 bin are counted and the fiducial color is 
applied according to their number densities (sec the anno- 
tated color wedge). The solid gray line in the lower panel 



http:/ /isc. astro. cornell.edu/~sloan/library/swsatlas/atlas. html 
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Fig. 8. The (Sll — L15) vs (N3 — Sll) color-color diagram of AKARI/IRC sources (upper panel) and ISO SWS sources (lower 
panel). In the AKARI panel only sources with S/N > 5 are included, and the (Sll — L15, N3 — Sll) plane is binned by 0.05 X 
0.05 mae 2 to calculate number of sources in each bin. The number densitv levels are logarithmic (sec the wcdecV The meanings 



12 



Y. Ita et al. 



shows the blackbody locus with the temperature ranging 
from 400 to 10000K in the IRC system. 

It is obvious that a wide variety of sources are detected 
in our survey. Comparing the two diagrams we can sta- 
tistically separate various types of objects. Two distinct 
main groups appear in the AKARI diagram: one centered 
around (Sll - L15, N3 - Sll) of (0.1, 1.0) and the other 
centered around (1.0, 4.0). 

The former group includes two peaks, one centered 
around (0.2, 0.2) and the other centered around (0.25, 
0.8). Based on a comparison with the ISO SWS clas- 
sification, the sources in the first peak are stars without 
circmstellar dust, and the second peak is associated to red 
giants with circmstellar dust emission. ISO SWS spectra 
indicate that red giants with O-rich dust as well as C-rich 
dust may once become blue in the (Sll — L15) color, and 
then turn red again in the course of their evolution. This 
can be explained by the change in strength of silicate or 
silicon carbide dust features in Sll band relative to the 
L15 flux. 

According to the ISO SWS spectra, the red group may 
contain various types of objects such as PPNe, PNc, 
OH/IR stars, and YSOs. Also, the background galaxies 
may have similar colors (see the following), which makes 
the classification more difficult. The NP spectroscopic 
data, and follow-up spectroscopic observations are needed 
to make clear classification of these objects. 

6.2. Color-magnitude diagrams 

An advantage to study sources in the LMC is that we 
can build color-magnitude diagrams (CMD) and estimate 
their absolute magnitudes by adding the distance modu- 
lus, which is reasonably well known (e.g., Feast & Walker 
1987). Based on the Spitzer SAGE catalog, Mcixncr et al. 
(2006) and Blum et al. (2006) presented infrared color- 
magnitude diagrams of LMC sources. Here we add Sll 
and L15 data, which are unique to the AKARI survey. 
We use intriguing combinations of the IRC, IRAC&MIPS 
and IRFS/SIRIUS bands and the results arc shown in 
Figure 9. Units of vertical axes are the apparent magni- 
tudes. It can be scaled to the absolute magnitude by sub- 
tracting the distance modulus. The corresponding wave- 
bands for the vertical axes are indicated at the top of 
each panel. The employed colors for the horizontal axes 
are indicated at the bottom of each panel. All of the 
color-magnitude planes are binned by 0.1 x 0.1 mag 2 to 
calculate the number of sources in each bin, and the fidu- 
cial color is given according to the number density levels in 
a logarithmic scale (see the wedge at the bottom). Other 
than the trends pointed out in Meixner et al. (2006) and 
Blum et al. (2006), new features are seen in our figure, 
primarily due to the addition of the Sll band data. 
6.2.1. New sequence for faint red giants with circumstel- 
lar dust 

Ramdani & Jorissen (2001) showed that the ratio of 
the ISO 11.5 fxm to DENIS K s 2 fxm flux density is a 
good indicator of dust mass-loss from red giants. The (K s 
— Sll) vs K s panel indicates which red giants show cir- 
cumstellar dust emission. It is clear that excess in Sll 



(K s — Sll > 0.5) is seen not only among the sources 
brighter than the tip of the first red giant branch (TRGB, 
K s ~ 12.5mag, Cioni et al. 2000b), but also among the 
sources below the TRGB. Since the brighter sources ex- 
ceed the K s luminosity of the TRGB, they should be an 
intermediate-age population and/or metal-rich old popu- 
lation AGB stars. On the other hand, the interpretation 
of the fainter sources is difficult. They can be metal-poor 
and old AGB stars that do not exceed the TRGB lumi- 
nosity, or red giants on the first red giant branch (RGB 
star). It is well known that the number ratio of AGB to 
RGB stars near the TRGB should be about 1/3 for the 
intermediate age stars (Renzini 1992). Thus, there are 
more RGB populations below the TRGB. The evidence 
for the mass-loss from the fainter sources is also seen in 
other panels. We see a bristle-shaped feature in (N3 — 
Sll) vs Sll panel, which is indicated by the arrow. The 
corresponding feature is marginal but also seen in the (S7 
— Sll) panel (indicated by the arrow). 

Lebzelter et al. (2006) obtained low-resolution mid- 
infrared (7.6 - 21.7 fxm) spectra of a star (V13) in the 
NGC104 with Spitzer. The K s band luminosity of the star 
is fainter than the TRGB luminosity of NGC104. They 
showed that the star is devoid of the 9.7 /xm emission 
band feature of amorphous silicate, but it has broad emis- 
sion features at 11.5 /urn (likely to be AI2O3) and 13 /zm 
(likely to be an Al— O stretching vibration). Aluminium 
oxide features have been detected from low mass-loss 
rate oxygen-rich red giants (Onaka et al. 1989; Kozasa 
& Sogawa 1997). Ita et al. (2007) suggested that the dust 
composition of V13 may be different from that of usual 
mass- losing AGB star (i.e., the brighter group). As it can 
be seen in Figure 2, the IRC Sll band includes all of these 
emission features. Therefore, the feature indicated by the 
arrow in the (N3 - Sll) vs Sll panel may be attributed 
to the red giants with aluminium oxide dust but without 
the silicate feature. Blum et al. (2006) identified a group 
of red but faint O-rich red giants in their ([8.0] — [24]) 
vs [8.0] color-magnitude diagram. The sources with Sll 
excess found in our (N3 — Sll) vs Sll panel may be sim- 
ilar to the red but faint O-rich red giants in Blum et al. 
(2006). Spectroscopic follow-up observations to identify 
the Sll excess with the AI2O3 band would be interesting. 
6.2.2. Carbon stars, dusty red giants, and YSO candi- 
dates on the CMDs 

We identify carbon stars, dusty red giants, and YSO 
candidates by the following criteria, and estimate their lo- 
cation on the CMDs. The results are shown in Figures 10 
and 11. 

• Optical carbon stars: As in section 5.4.3, we cross- 
identified our catalog with the carbon star catalog 
in the LMC (Kontizas et al. 2001) based on optical 
prism spectroscopy. They are shown by the green 
dots in Figures 10 and 11. 

• Dusty red giants: We select sources that satisfy the 
same conditions employed in section 5.4.3. We use 
the term "dusty red giants" here, but a significant 
fraction of the sources in this category can be ac- 




Fig. 9. Color magnitude diagrams of AKARI LMC sources using several combinations of IRC, IRAC, and IRSF/SIRIUS bands. 
The vertical axis is in apparent magnitude at the corresponding wavebands, which arc indicated at the top of each panel. The allows 
indicate the newly found feature, which can be attributed to the red giants that have luminosities below the tip of the first red giant 
branch. 



Table 3. Assumed zero magnitude fluxes and their corresponding reference wavelengths for the existing catalogs. 



Wavebands 


Zero magnitude flux 

[Jy] 


Reference wavelength 

H 


Reference 


U 


1649 


0.3745 


Cohen et al. (2003a) 


B 


4060 


0.4481 




V 


3723 


0.5423 


55 


I 


2459 


0.8071 


J? 


J 


1594 


1.235 


Cohen, Wheaton, & Megeath (2003b) 


H 


1024 


1.662 




K s 


666.8 


2.159 


» 


[3.6] 


280.9 


3.550 


IRAC Data Handbook (2006) 


[4.5] 


179.7 


4.493 




[5.8] 


115.0 


5.731 


?? 


[8.0] 


64.1 


7.872 


?? 


[24] 


7.14 


23.68 


MIPS Data Handbook (2007) 
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Fig. 10. The same as Figure 9, but with spectroscopically confirmed carbon stars (green dots) by optical survey, candidates of 
dusty red giants (red dots), and candidates of YSOs (blue dots) are identified. 



tually infrared carbon stars, which are elusive in 
the above optical spectroscopic survey. According 
to Nikolacv & Weinberg (2000), stars with (J — K s ) 

> 2.0 in the LMC are primarily carbon-rich AGB 
stars. This can be explained as the J band flux of 
carbon-rich mass-losing red giants are attenuated by 
C2 and CN absorptions (Loidl et al. 2001). This cri- 
terion may be satisfied also by some dusty oxygen- 
rich AGB stars, but their number should be small 
compared to carbon stars. To confirm their circum- 
stellar chemistry, our NP spectroscopic data would 
be effective. Sources in this category are shown by 
the red dots in Figures 10 and 11. 

• YSO candidates: Based on the color-color diagram 
(Figure 8), we select sources that satisfy (N3 — Sll) 

> 2.0 and (Sll - L15) > 1.0. As seen in the ISO 
SWS panel of Figure 8, sources selected by this cri- 
teria include not only YSOs but also other popu- 
lations such as PPNe and PNe. Also, background 
galaxies may have colors similar to sources in this 
category (see the following). Sources in this cate- 
gory are shown by the blue dots in Figures 10 and 
11. 



Also, we calculated bolometric magnitudes for sources 
that satisfy the following conditions. 

• Sources detected at least in two of four optical (U, 
B, V, and I) bands. 

• Sources detected at least in three of seven near- 
infrared (J, H, K s , N3, [3.6], [4.5], and [5.8]) bands. 

• Sources detected at least in three of four mid- 
infrared (S7, [8.0], Sll, and L15) bands. 

• Sources detected at L24 and/or [24] bands. 

Table 3 summarizes the assumed zero magnitude fluxes 
and their corresponding reference wavelengths that we 
used to calculate the flux densities. For the IRC bands, we 
use the zero magnitude fluxes and reference wavelengths 
listed in Table 1. Since the zero magnitude fluxes for 
the IRSF/SIRIUS JHK S bands are not officially deter- 
mined, we converted the IRSF system magnitudes into 
the 2MASS system by using the conversion equations pro- 
vided in Kato et al. (2007) before calculating the flux 
densities. Then we use a cubic spline to interpolate the 
spectral energy distributions and integrate them from the 
shortest wavelength at which the flux is available to 24 fan 
to obtain the apparent bolometric magnitudes. Finally, 
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Fig. 11. The same as Figure 10 but with the vertical axis in absolute bolometric magnitude. We assumed a distance modulus of 
18.5 mag to the Large Magellanic Cloud. The meanings of colors of marks are the same in Figure 10. 



we scale them to the absolute magnitudes by assuming 
a distance modulus of 18.5 mag to the LMC. Figure 11 
shows the result. For very red sources, such as galaxies, 
planetary nebulae, and extremely dusty AGB stars, the 
bolometric magnitudes may be underestimated to a large 
extent because the fluxes longward of 24 /im are not in- 
cluded. Therefore the bolometric magnitudes discussed 
here should be only lower limits. 

It is apparent that optical carbon stars, dusty red gi- 
ants, and YSO candidates are well separated from each 
other by color and luminosity. Except for the ncar- 
infrarcd CMD, YSO candidates have redder colors than 
that of the others, and dusty red giants are brighter in 
the mid- infrared bands than optical carbon stars. In the 
(N3 — [5.8]) vs [5.8] CMD, we see a group of sources 
that distribute along the vertical sequence around (N3 
— [5.8] ~ 0.0). These sources are likely to be foreground 
Galactic stars that are classified as dwarf or giant. There 
is another feature in the diagram. We see a sequence for 
dusty red giants (red dots), which grows toward the red- 



der and brighter part of the CMD. Above the sequence, 
we see another branch that grows toward the reddder 
and brighter part of the CMD. The brighter branch is 
clear in (N3 — [5.8]) vs [5.8] diagram, and is also seen 
in the (N3 — Sll) vs Sll diagram. The corresponding 
bolometric magnitude CMDs tell us that most of dusty 
red giants have absolute bolometric magnitudes (M^ i) 
ranging from —4 to —5, while the sources on the brighter 
branch have Mi^i brighter than —6 mag and they are sat- 
urated in the IRSF/SIRIUS bands. It is suggested that 
sources on the brighter branch have stellar masses heav- 
ier than those of the sources on the sequence for dusty 
red giants. Their bolometric magnitudes agree well with 
the bolometric magnitudes for 88 red supergiants in the 
LMC (Oestreicher & Schmidt-Kaler 1998). Therefore we 
conclude that the sources on the brighter branch are red 
supergiants. 

As we suggested in section 6.1.2, the (Sll — LI 5) color 
of dusty red giants once becomes blue during their evolu- 
tion due to the growth of silicate or silicon carbide emis- 
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sion in the Sll band. In fact most dusty red giants (red 
dots), if not all, are bluer in (Sll — L15) than optical 
carbon stars (green dots). 

Table 4. Galactic coordinates of reference fields. 



ID 


latitude 


longitude 




[degree] 


[degree] 


1 


86.26824 


34.38021 


2 


87.56487 


34.71097 


3 


87.67876 


35.02329 


4 


244.41766 


-34.71058 


5 


68.25446 


35.19998 
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Fig. 12. Color-magnitude diagrams for sources in reference 
fields. 

6.2.3. Background galaxies and foreground Milky Way 
sources 

To estimate how many background galaxies and fore- 
ground sources are expected in the CMDs, we use the 
leftover NIR and MIRS data that were kindly offered by 
the AKARI stellar working group (PI. Y.Nakada). The 
data were obtained when they took MIRL images of glob- 
ular clusters. MIRL observations concurrently yield NIR 
and MIRS images that observe 25' away from the MIRL 
center (= target center). Also they observed some local 
group galaxies and took several reference images for them, 
which were about 1 degree away from the target centers. 
These reference data were also included in the following 
analysis. Among those offered by AKARI stellar working 
group, we select fields whose absolute value of galactic 
latitudes \b\ are in the range of 30° < |6| < 40° that are 
comparable to the mean value of the AKARI LMC sur- 
vey regions (\b\ ~ 32.5°). We reduce the data in the same 
way as we did for the LMC survey. The total area of the 
reference fields is 500 arcmin 2 . The coordinates of the ref- 
erenced fields are listed in Table 4. We check IRAS 24, 60 
and 100 ^m images of the corresponding coordinates, and 
found no appreciable indication of star forming regions 
and molecular clouds. 



Using the data of the reference fields, we make CMDs 
that can be directly compared to Figures 9 and 10. They 
are shown in Figure 12. Most, if not all, of the sources 
in Figure 12 arc expected to be background galaxies and 
foreground Milky Way stars, with a few possibility of stars 
in the outer skirts of target globular clusters. We check 
the tidal radii of the target globular clusters (Harris 1996), 
and they are about a half of the offset value of 25'. In the 
following analysis, we will ignore the small contribution 
from the stars in the outer skirts of each globular cluster. 
Due to the high source density, in particular in the cen- 
tral part of the LMC, the completeness limit of the LMC 
survey data must be shallower than that of the reference 
fields. Therefore, we only count the number of sources 
that satisfy the following conditions. 

(1) For presumable background galaxy: 

• (N3 - Sll) > 2 [mag] 

• 10 < Sll [mag] < 12.2 

(2) For presumable foreground stars: 



(N3 - Sll) < 2 [mag] 
9 < Sll [mag] < 12.2 



Then we found 72 background galaxies and 26 fore- 
ground stars in the reference fields. The differences in 
number counts for both populations among the reference 
fields are small. These results imply that we can expect 
about 500 background galaxies and 190 foreground stars 
in a square degree. This estimation agrees with the Spitzer 
results (Fazio et al. 2004b, Meixner et al. 2006). On the 
other hand, we found 4483 and 14951 sources in our cat- 
alog (compiled from 10 deg 2 survey data) that satisfy the 
above conditions, respectively. Therefore most, if not all, 
of the sources that satisfy the condition (1) in our catalog 
should be background galaxies. Meanwhile, the fraction 
of foreground stars in the sources that satisfy condition 
(2) is estimated to be about 13%. 

7. Summary 

We carried out an imaging and spectroscopic survey of 
a 10 deg 2 area of the LMC using IRC onboard AKARI. 
In this paper we discuss the imaging data obtained in 
the survey and describe a preliminary photometric cata- 
log of bright point sources. The first point source cata- 
log is planned to be released to the public in 2009. We 
cross-identified our catalog with the existing optical, near- 
infrared, and mid-infrared photometry catalogs, and used 
the combined data to build color-magnitude diagrams us- 
ing several combinations of interesting wavebands and also 
to calculate absolute bolometric magnitudes. By virtue 
of the Sll and L15 bands that are unique to IRC, we 
found new interesting features in the color-magnitude di- 
agrams. The present IRC data and the upcoming results 
from AKARI All-Sky survey together with existing radio 
and near-infrared ground-based observational results will 
provide a significant database to study the star-formation 
history and material circulation within a galaxy. 
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